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ANALYSIS  OF  LINEAR  PARAMETER  VARYING  SYSTEM  MODELS  BASED  ON 

REACHABLE  SETS 


JONG-YEOB  SHIN  * 

Abstract.  This  paper  presents  the  analysis  method  of  quasi-LPV  models,  comparing  the  ellipsoid  set 
which  contains  the  reachable  set  of  a  nonlinear  system  to  define  which  quasi-LPV  model  is  less  conservative 
to  represent  the  nonlinear  dynamics.  Three  quasi-LPV  models  are  constructed  from  a  nonlinear  model 
using  three  different  methods,  to  facilitate  synthesis  of  an  LPV  controller  for  the  nonlinear  system.  The 
comparison  results  of  closed-loop  system  performance  with  synthesized  LPV  controllers  correspond  to  the 
analysis  results  of  quasi-LPV  models. 

Key  words.  LPV  model,  reachable  set,  LPV  control  synthesis 

Subject  classification.  Guidance  and  Control 

1.  Introduction.  One  of  control  schemes  for  nonlinear  systems  is  a  linear  parameter  varying  (LPV) 
technique  [7,  6,  11,  4,  3,  10].  This  approach  is  particularly  appealing  in  that  nonlinear  plants  are  treated 
as  linear  systems  with  varying  parameters.  This  allows  linear  control  techniques  to  be  applied  to  nonlinear 
systems.  In  applying  LPV  analysis  and  synthesis  methodology  to  nonlinear  systems,  an  LPV  plant  model 
of  a  nonlinear  system  is  required  to  describe  the  nonlinear  dynamics.  Since  an  LPV  controller  is  synthesized 
based  on  an  LPV  model  and  is  applied  to  control  the  nonlinear  system,  it  is  important  to  choose  which  LPV 
model  is  used  for  an  LPV  controller  synthesis  to  lead  a  less  conservative  result. 

There  are  three  different  approaches  to  generate  an  LPV  model  from  a  nonlinear  mathematical  model 
of  a  nonlinear  system.  Conventionally,  an  LPV  model  is  constructed  by  the  set  of  linearized  models  around 
equilibrium  points.  Another  approach  is  state  transformation  which  changes  state  coordinates  to  remove 
nonlinearity  in  the  dynamics  [5,  8].  The  other  approach  is  function  substitution  which  formulates  nonlinear 
functions  into  quasi-LPV  form  functions  [15,  16,  9].  The  LPV  models  generated  by  these  three  methods  have 
been  discussed  in  terms  of  accuracy  to  present  a  nonlinear  dynamics  by  comparing  time  simulation  results 
with  the  pre-defined  input  signals  [8,  1].  Comparison  of  the  time  responses  is  one  of  approaches  to  decide 
which  LPV  model  will  be  used  for  LPV  analysis  and  synthesis  of  a  nonlinear  system. 

There  are  possibilities  for  existence  for  different  LPV  models  to  produce  time  responses  of  a  nonlinear 
system  accurately  within  ignorable  error  range.  Therefore,  in  this  case,  comparison  of  the  time  simulation 
results  of  the  LPV  models  is  not  sufficient  to  analyze  the  models.  In  this  paper,  one  of  approaches  is 
demonstrated  to  analyze  LPV  models  to  find  which  LPV  model  is  less  conservative  to  describe  the  nonlinear 
system. 

In  this  paper,  three  different  LPV  models  of  a  nonlinear  system  are  provided  by  three  methods  and  time 
simulation  of  three  LPV  models  with  the  pre-defined  input  signals  are  presented.  Ellipsoids  which  contain 
the  reachable  set  of  the  nonlinear  system  are  calculated  according  to  each  LPV  model  and  compared  the 
sizes  of  the  ellipsoids.  Also,  each  LPV  controller  is  synthesized,  based  on  each  LPV  model  and  is  applied  to 
control  the  nonlinear  system. 

*ICASE,  MS.  132c,  NASA  Langley  Research  Center,  Hampton,  VA,  23681  (emaihj. y.shin@larc.nasa.gov).  This  research  was 
supported  by  the  National  Aeronautics  and  Space  Administration  under  NASA  Contract  No.  NASl-97046  while  the  author 
was  in  residence  at  ICASE,  NASA  Langley  Research  Center,  Hampton,  VA  23681-0001 
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Outline  of  this  paper  is  follows.  In  section  2,  conventional  LPV  control  synthesis  used  in  this  paper  is 
briefly  summarized.  In  section  3,  an  analysis  methodology  of  quasi-LPV  models  is  presented,  calculating 
an  ellipsoid  which  contains  a  reachable  set.  In  section  4,  an  example  of  analysis  of  quasi-LPV  models  is 
demonstrated  and  also  nonlinear  simulations  with  LPV  controllers  synthesized  based  on  each  quasi-LPV 
model  are  presented  for  comparison.  This  paper  concludes  with  a  brief  summary  in  section  5. 

2.  LPV  Control  Synthesis.  In  this  section,  a  quasi-LPV  system  is  deflned  and  an  LPV  control 
synthesis  methodology  is  briefly  described.  Consider  a  generalized  linear  open-loop  system  as  functions 
of  parameters  p(t)  G  V.  For  a  compact  subset  V  C  71^ ,  the  parameter  variation  set  denotes  the  set  of 
all  piecewise  continuous  functions  mapping  R  (time)  into  V  with  a  flnite  number  of  discontinuities  in  any 
interval,  where  s  is  number  of  parameters.  An  LPV  open-loop  system  can  be  written  as 


x{t) 

'A{p{t))  Bi{p{t))  B2{p{t))' 

x{t) 

e{t) 

= 

Ci{p{t))  0  Di2(p(i)) 

d{t) 

C2{p{t))  D2l{p{t))  0  . 

u{t)_ 

where  y{t),  e(t),  d{t)  and  u{t)  are  measurements,  errors,  disturbances,  and  control  signals.  A  quasi-LPV 
system  is  deflned  when  scheduling  parameter  vector  p  contains  part  of  state  vector  as  [xi  pe],  where  pe  is 
an  exogenous  scheduling  parameter  vector  and  the  state  vector  x  =  [xi  X2].  Hereafter,  p  denotes  p{t).  The 
induced  C2  norm  of  d  to  e  is  deflned  as 


sup 


Mk 

\\d\W 


In  an  LPV  synthesis  methodology,  suppose  there  is  an  LPV  output  feedback  controller  K{p)  which 
stabilizes  the  closed-loop  system  exponentially  and  makes  the  induced  >C2-norm  of  d  to  e  less  than  7.  The 
controller  K{p)  can  be  written  as 

Xk 
u 


M  {p) 
Ckip) 


Skip) 

Dk  {p) 


Xk 

_y  _ 

(2.2) 


An  LPV  controller  K{p)  can  be  constructed  from  solutions  of  V  G  and  Y  G  of  the  following 

optimization  problem  [4]. 


subject  to 


min  7 

x,Yen^^^ 


(2.3) 


'xA^{p)+A{p)X-B2{p)Bj{p) 

xc^iip) 

l-^B,{p) 

Cii{p)X 

0 

<  0 

(2.4) 

0 

~dnd 

'A{p)Y  +  YA^{p)-C^{p)C2{p) 

yBii{p) 

{p) 

BiAp)y 

0 

<0, 

(2.5) 

7~^Ci{p) 

0 

-In. 
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X  >  0,  y  >  0 


where 


A(p)  =  A(p)  -  B2(p)Ci2ip),  A{p)  =  A{p)  -  Bi2ip)C2ip),  (2.7) 

and  n  is  number  of  states  of  the  generalized  open-loop  system.  Note  that  X  and  Y  are  constant  positive 
definite  matrices. 

A  method  to  construct  an  LPV  controller  K{p)  from  the  solution  matrix  X  and  Y  of  the  LMI  optimiza¬ 


tion  problem  is  taken  from  Ref.  [4].  An  LPV  controller  is  constructed  as  [4]: 

Ak(p)  =  A{p)  +  B2{p)F{p)  +  Q-^YLip)C2ip)  -  j-^Q-^M{p),  (2.8) 

Bkip)  =  -Q-^YL{p),  (2.9) 

Ck{p)=F{p),  (2.10) 

Dk{p)=0  (2.11) 


where  matrices  Q,  F(p),  L(p),  and  M{p)  are  defined  as 

Q  =  Y-j-^X-\ 

F{p)  =  -[B^{p)X-^+Df,{p)C^{p)], 

L{p)  =  -[Y-^C^{p)+B,{p)D^,{p)], 

M{p)=H{p)+^^Q[-Q-^YL{p)D2i{p)-Bi{p)]Bl{p)X-\ 

Matrix  H  {p)  is  defined  as 

H{p)  =  -[X-^AF{p)+AF{pfX-^+Cl{p)CF{p)+l-^X-^Bi{p)Bi{pfX-^] 

with  Af{p)  =  A{p)  -h  B2{p)F{p)  and  Cf{p)  =  -h  Di2F{p).  The  closed-loop  system  with  the  controller 
K{p)  is  exponentially  stable  and  the  induced  C2  norm  is  less  than  7.  The  proof  can  be  found  in  Ref.  [4]. 
To  make  the  optimization  problem  of  equation  (2.3)  computationally  tractable,  scheduling  parameters  p  are 
discretized  into  grid  points.  Thus,  infinite  LMI  constraints,  equations  (2.4)-(2.5),  are  presented  as  finite 
number  of  LMI  constraints.  Note  that  the  LPV  controller  may  be  different  based  on  scheduling-parameter 
grid  points. 

Note  that  the  open-loop  LPV  system  matrices  A(p),  B{p)^  C{p)  and  D{p)  are  used  to  construct  the 
LPV  controller  (see  equations  (2.8)-(2.11)).  The  LPV  controller  can  be  different  based  on  the  different  LPV 
models  of  a  nonlinear  system,  even  if  the  different  LPV  models  can  produce  the  exact  same  input-output 
time  responses  for  each  other. 

3.  Analysis  of  Quasi-LPV  Model.  In  this  section,  one  of  approaches  to  analyze  an  LPV  model  is 
presented  in  terms  of  the  sizes  of  ellipsoids  which  contain  a  reachable  set  of  a  nonlinear  system.  Before  we 
introduce  an  analysis  method  of  LPV  models  of  a  nonlinear  system,  a  class  of  nonlinear  systems  used  in  this 
paper  is  defined  as  follows. 

Consider  a  nonlinear  system  in  which  an  input  vector  enters  affinely.  A  nonlinear  system  can  be  written 
as 


X  =  F{x)  -h  B{x)u 


(3.1) 
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where  a  state  vector  x  is  in  7^^,  an  input  vector  u  is  in  7^"^,  continuous  function  F{x)  :  71^  and 

B{x)  :  VF'  reachable  set  of  a  nonlinear  system  with  bounded-energy  inputs  is  defined  by 


= 


x{T) 


X  =  F{x)  -h  B{x)u,  x(0)  =  0, 


u^udt  < 


T  >  0 


(3.2) 


where  a  is  a  given  positive  constant.  The  reachable  set  of  a  nonlinear  system  is  bounded.  Here  an  initial 
point  of  X  is  defined  as  0  without  loss  of  generality.  Definition  of  an  invariant  set  used  in  this  paper  is  follows: 


Definition  3.1  Invariant  set 

Let  a  set  e  denote  a  set  eentered  at  the  origin 

s  =  {xen^\v{x)  <  1}. 

The  set  e  is  said  to  he  invariant  if  for  every  trajeetory  x{t)  of  a  nonlinear  system,  x(0)  G  e  implies  x{t)  G  e 
for  all  t.  A  funetion  V  (x)  is  a  Lyapunov  funetion. 

It  is  easily  shown  that  when  V{x)  >0  and  <  0  for  all  trajectory  satisfying  the  nonlinear  dynamics  in 

equation  (3.1),  the  set  e  is  invariant. 

A  quasi-LPV  model  of  a  nonlinear  system  can  be  produced  with  or  without  a  bounded  uncertainty  block 
to  capture  nonlinear  dynamics. 

•  Case  1  :  a  quasi-LPV  model  can  describe  the  nonlinear  dynamic  model  without  an  uncertainty 
block.  A  quasi-LPV  model  is 


X  =  A{x)x  -h  B{x)u  =  F{x)  -h  B{x)u. 


(3.3) 


The  reachable  set  of  a  quasi-LPV  model  is  defined  as 


Fbe  = 


X 


X  =  A{x)x  -h  B{x)u,  x(0)  =  0, 


u^udt  <  , 


T  >  0 


(3.4) 


Thus,  it  is  obvious  that  the  reachable  set  IZbe  is  equal  to  IZni- 
•  Case  2:  a  quasi-LPV  model  can  describe  the  nonlinear  dynamic  model  with  a  bounded  uncertainty 
block.  A  quasi-LPV  model  is 


X  =  A{x)x  -h  Bu{x)u  -h  Bw{x)w, 
z  —  C 2,{x^x  “h  D 2,uix^u  “h  D ^.wix^uj , 
w  =  Az 


(3.5) 

(3.6) 

(3.7) 


where  B^{x)  e  TZ^^p,  w  e  TZp,  C,  e  e  e  TZ^^p,  and  ||A||  <  /?.  The 

reachable  set  of  a  quasi-LPV  model  is  defined  as 


IZb 


eu 


<  X 


X  =  A{x)x -\- B{x)u -\- B^{x)w, 
z  =  Cz{x)x  -h  Dzu{x)u  -h  Dzw{x)w, 
w  =  Az 


^(0)  =  0,  f 

Jo 


mudt  <  ,  T  >  0 


>  . 


(3.8) 


Assume  that  there  exists  A  such  that 


hZni 


F  IZbeu- 


(3.9) 
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Note  that  calculating  the  size  and  structure  of  A  to  validate  equation  (3.9)  is  out  of  this  paper  scope. 
Suppose  there  exists  a  Lyapunov  function  V  such  that 


dV(x(t)) 

dt 


(3.10) 


for  every  x{t)  and  u{t)  satisfying  equation  (3.3)  or  equations  (3.5)-(3.7).  Then  there  exists  the  invariant  set 
{x  I  V{x)  <  1}  which  contains  the  reachable  set  TZte  or  TZteu,  according  to  each  quasi-LPV  model. 


3.1.  Singular  Quadratic  Lyapunov  Function.  In  this  section,  an  LMI  optimization  problem  is 
formulated  to  calculate  the  smallest  ellipsoid  which  contains  the  reachable  set,  using  a  singular  quadratic 
Lyapunov  function  V{x)  =  x^Px.  For  the  two  cases:  Case  1  and  Case  2,  the  LMI  optimization  problems 
are  formulated  as  follows. 

For  case  1: 


sup  Tr{P)  (3.11) 

p^-JlnXn 


subject  to 


A{x)^P  +  PA{x)  PB{x) 

BixfP  -^I 


<  0. 


(3.12) 


The  LMI  constraint  of  equation  (3.12)  is  easily  derived  from  equation  (3.3).  Thus, 

the  set  {x{t)\x{t)^ Px{t)  <  l,t  >  T}  contains  the  reachable  set  Tl^e.  Also,  it  is  noticed  from  the  LMI 
constraint  that  the  ellipsoid  is  an  invariant  set. 

For  case  2: 


suprr(P)  (3.13) 

p 


A^{x)P  +  PA{x)  +  tCj  {x)Cz{x) 
Bl{x)PptDU^)Cz{x) 

.  BUx)P  +  tDUx)Cz{x) 


PBu{x)ptCj{x)D 

zu  (x) 

+  tD"^^(yX)D  zu{x) 


0 


PBn,{x)  +  tCj  {x)Dzn,{x) 
0 

zw{x)  —  -^I 


<  0 


(3.14) 


P  >  0,  t>0  (3.15) 

where  ||A||  <  [3  and  [3  is  given  constant.  The  LMI  constraint  of  equation  (3.14)  is  easily  derived  from 
—  ■^Il'^(^)ll2  and  equations  (3.5)-(3.7),  using  5-procedure  [14].  The  set  {x{t)\x{t)'^  Px{t)  <  l,t  >  T} 
contain  the  reachable  set  IZheu-  Also,  it  is  noticed  from  the  LMI  constraint  that  the  ellipsoid  is  an  invariant 
set. 

The  size  of  the  ellipsoid  is  defined  as 


n 

5p  =  CnV(P)  (3.16) 

i 

where  C  is  constant  which  is  dependent  on  the  geometry  of  the  ellipsoid.  Thus,  it  is  easy  to  compare  the 
size  of  the  ellipsoid  set  {x\x^ Px  <  1}. 
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3.2.  Parameter  Dependent  Lyapunov  Function.  In  this  section,  the  LMI  optimization  problem  is 
formulated  with  the  parameter-dependent  Lyapunov  function  V (x)  =  x^P(x)x.  LMI  constraints  are  written 
for  each  case  in  the  same  manner  of  formulation  of  LMI  constraints  in  equations  (3.12)  and  (3.14). 

For  case  1: 


subject  to 


For  case  2: 


sup  Tr{P{x)) 


'A{x)'^P{x)  +  PA{x)  -h  P{x)  P{x)B{x)) 
B{xYP{x)  -^I 


<  0. 


(3.17) 


(3.18) 


sup  Tr{P{x)) 

P(a?):7^^^7^^xn 


(3.19) 


A^{x)P{x)  +  P{x)A{x)  +  tc^  {x)Cz{x)  +  P{x)  P{x)Bu{x)  +  tC^  {x)Dzu{x) 
Bu  {x)P{x)  +  {x)Cz  {x)  “  {x)Dzu  {x) 

Bl(x)P(x)  PtD'^^(x)Cz(x)  0 


P{x)Byj  (x)  +  tCj  {x)Dzw  {x) 

0 

tD^^{x)Dzw{x)  — 


<  0  (3.20) 


P  >  0,  t  >  0. 


(3.21) 


Using  the  LMI  optimization,  the  smallest  ellipsoid  to  capture  the  reachable  set  of  the  nonlinear  dynamics 
can  be  calculated.  To  solve  the  LMI  optimization  problem,  the  basis  functions  of  P{x)  are  required.  The 
details  of  solving  the  LMI  with  P{x)  will  be  explained  in  the  next  section.  Note  that  the  size  of  the  set 
{x\x{t)'^ P{x{t))x{t)  <  1}  is  not  easy  to  calculate  since  it  is  not  ellipsoid. 


4.  Example.  In  this  section,  LPV  models  of  a  nonlinear  system  are  generated  by  three  different  meth¬ 
ods  (Jacobian  linearization,  state  transformation,  and  function  substitution)  and  each  LPV  model  is  simu¬ 
lated  with  pre-defined  input  signals.  The  simulation  results  are  compared  with  nonlinear  simulation  results 
to  notice  that  quasi-LPV  models  can  capture  the  nonlinear  dynamics.  Also,  the  size  of  an  ellipsoid  which 
contains  the  reachable  set  of  the  nonlinear  system  is  calculated  for  each  LPV  model.  An  LPV  controller  is 
synthesized  based  on  each  LPV  model,  respectively  and  is  simulated  with  nonlinear  system  to  compare  the 
closed-loop  performance. 

A  nonlinear  system  taken  from  Ref.  [2]  is 


Xi 

X2 


-Xi 

Xi  -  \x2\x2  -  10 


+ 


1 

0 


u,  y  =  X2. 


(4.1) 


It  is  noted  that  systems  with  similar  types  of  nonlinearity  are  frequently  encounted  in  practice  [2,  13]. 


4.1.  Quasi-LPV  Models.  Using  Jacobian  linearization  around  trim  points,  the  set  of  linearized  models 
can  present  a  quasi-  LPV  model  of  the  nonlinear  system.  The  linearized  model  at  a  trim  point  is 


5xi 

5X2 

Sxi 

Sx2 


+ 


1 

0 


Su,  Sy  =  (5x2, 


(4.2) 


where  Sxi  =  xi  —  xiq,  6x2  =  X2  —  X2o,  Su  =  u  —  uq,  and  Sy  =  y  —  yo.  The  associated  quasi-LPV  system  is 


m 


-1 

1 


0 

-2z 


m 


1 


+ 


0 


u, 


m 


(4.3) 
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where  the  range  of  2:  =  |x2o|  is  defined  arbitrary  from  0  to  5.  It  is  noted  that  states,  r]i  and  7^2,  are  defined 
as  deviation  from  each  trim  point  {xiq,X2q,uo).  Thus,  when  scheduling  parameter  2:  changes  the  definitions 
of  states  of  the  LPV  model  are  changed. 

To  use  function  substitution  [15,  16,  9],  the  nonlinear  system  is  rewritten  as 


where 


Xl 

’-1  o' 

Xl 

+ 

'i' 

—  XIq  -h  Uq 

= 

u  -h 

_X2_ 

1  0 

_X2_ 

0 

Xi„  -  \X2  +  X2o  1  {X2  +X2o)  -  10_ 

(4.4) 


Xi  =  Xi  -  Xio  ,  X2  =  X2  -  X2o  ,  U  =  U  -  Uq.  (4.5) 

The  nonlinearity  in  equation  (4.4)  is  substituted  for  a  function  in  quasi-LPV  form. 


/(^2)  = 

A  quasi-LPV  model  is  written  as 


[1^20  k2o  -  1^2  +^2ol(^2  +^2o)]/^2,  ^2  ^  0, 

0,  X2  =  0. 


-1  0 
1  f{x2) 


+ 


(4.6) 


(4.7) 


Here,  a  trim  point  is  set  as  (xiq,  X2q)  =  (11,  1).  Note  that  the  definition  of  states  of  this  quasi-LPV  model 
is  deviation  from  the  trim  point.  Also,  the  definition  of  states  does  not  change  as  a  scheduling  parameter  X2 
changes. 

A  quasi-LPV  model  of  the  nonlinear  system  can  be  generated  by  changing  state  coordinates  [5] 


T]l=  xi-  xio(x2),  (4.8) 

m  =  X2  (4.9) 


where 


Xlo{x2)  =  \X2\X2  +  10. 


(4.10) 


A  quasi-LPV  model  is 


_ 

■-1-2|7?2|  o' 

+ 

'i' 

j]2_ 

1  0 

0 

(4.11) 


where  r  =  u  —  uq{x2).  Note  that  definition  of  control  input  r  and  state  r]i  are  changed  as  scheduling 
parameters  vary.  When  the  quasi-LPV  model  is  simulated,  the  variations  of  uo{x2)  and  771  should  be 
compensated  as  the  scheduling  parameter  X2  changes. 


4.2.  Quasi-LPV  Model  Simulations.  There  are  three  quasi-LPV  models  to  describe  the  nonlinear 
system.  To  compare  the  simulation  results  between  the  nonlinear  system  and  the  quasi-LPV  model,  a 
performance  index  J  is  introduced  as  C2  norm  of  error  in  finite  time  T, 

J  =  [  eitfeit)  dt=  f  {yni{t)  -  yLPv{t)Y{yni{t)  -  VLPvit))  dt  (4.12) 

Jo  Jo 

where  yni  and  yLPV  are  measurements  of  the  nonlinear  system  and  the  quasi-LPV  model,  respectively.  Here, 
T  is  set  as  30  sec  in  this  example. 
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The  quasi-LPV  models  are  simulated  for  the  two  different  input  signal  sets  (see  Figure  5.1),  respectively. 
The  simulation  results  in  Figure  5.2  show  that  the  quasi-LPV  models  present  the  nonlinear  dynamics  very 
accurately.  The  differences  of  time  responses  between  the  quasi-LPV  models  and  the  nonlinear  system  are 
not  noticeable  in  Figure  5.2.  The  performance  index,  J,  for  all  quasi-LPV  models  is  less  than  10“^  for 
both  different  input  cases.  The  indexes,  NL,  Q-LPVj,  Q-LPVg,  and  Q-LPV/,  in  Figure  5.2  are  denoted 
as  the  nonlinear  system,  the  quasi-LPV  model  by  Jacobian  linearization,  the  quasi-LPV  model  by  state- 
transformation,  and  the  quasi-LPV  model  by  function  substitution.  Note  that  all  quasi-LPV  models  can 
present  all  state  dynamics  of  the  nonlinear  system  very  accurately.  Therefore,  comparing  performance  index 
is  not  sufficient  to  choose  which  quasi-LPV  model  will  be  used  for  LPV  controller  synthesis. 

4.3.  Reachable  Set.  In  this  section,  the  smallest  ellipsoid  which  contains  the  reachable  set  of  the 
nonlinear  model  is  calculated  for  each  quasi-LPV  model. 

4.3.1.  Singular  Quadratic  Lyapunov  Function.  Assume  a  Lyapunov  function  V (x)  =  x^Px, 
P  G  is  used  to  calculate  the  smallest  set  {x\x^ Px  <  1}  which  contains  the  reachable  set  of  the 

nonlinear  system.  Consider  a  quasi-LPV  model  without  an  uncertainty  block,  which  can  capture  the  nonlin¬ 
ear  dynamics  (Case  1).  To  make  the  LMI  optimization  problem  in  equation  (3.11)  computationally  tractable, 
the  infinite  number  of  LMI  constraints  are  converted  into  the  finite  number  of  LMI  constraints  defined  each 
grid  points  over  the  parameter  space.  The  grid  points  in  this  example  are  defined  as 


X2,  G  {0.1, 0.5, 0.9,2, 3, 4,5}. 


(4.13) 


Also,  bounded  energy  inputs  are  required  to  solve  the  LMI  optimization  problems.  Hereafter,  the  bounded 
energy  inputs  are  defined  u^u  <  0.2  in  this  example. 

The  solution  matrix  P  of  the  LMI  optimization  problem  is  calculated,  based  on  the  quasi-LPV  models 
generated  by  Jacobian  linearization  and  function  substitution,  respectively.  The  ellipsoids  corresponding  to 
the  solution  matrix  P  are  shown  in  Figure  5.3  for  each  quasi-LPV  model.  In  Figure  5.3,  possible  trajectories 
of  the  nonlinear  system  with  pre-defined  several  input  signals  bounded  by  constant  energy  are  also  plotted 
as  dashed  lines.  It  is  observable  that  the  two  ellipsoids  capture  the  candidate  trajectories  to  represent  the 
reachable  set  of  the  nonlinear  system.  Also,  it  is  noticeable  that  the  size  of  the  ellipsoid  based  on  the  function 
substitution  quasi-LPV  model  (Q-LPV/)  is  smaller  than  that  of  the  Jacobian  linearization  quasi-LPV  model 
(Q-LPVj).  Thus,  the  quasi-LPV  model  generated  by  function  substitution  is  less  conservative  than  the 
quasi-LPV  model  generated  by  Jacobian  linearization  to  contain  the  reachable  set  of  the  nonlinear  system. 

Since  one  of  eigenvalues  of  the  quasi-LPV  model  generated  by  state  transformation  is  0  for  all  grid 
points,  it  is  not  computationally  tractable  to  calculate  the  size  of  the  ellipsoid  to  contain  the  reachable  set. 
It  is  noticeable  that  the  quasi-LPV  model  can  not  satisfy  one  of  the  conventional  assumptions  for  Poo  and 
LPV  synthesis  [4]  : 


a  bl 
c2  d21 


have  full  row  rank. 


(4.14) 


Thus,  this  analysis  method  can  not  apply  for  a  quasi-LPV  model  which  has  zero  eigenvalue  over  all  grid  points. 
Hereafter,  we  discuss  two  quasi-LPV  models  generated  by  Jacobian  linearization  and  function  substitution 
methods. 

Consider  that  a  quasi-LPV  model  with  an  uncertainty  block  to  represent  the  nonlinear  dynamics  (Case  2). 
Here,  unmodeled  dynamics  is  assumed  as  an  input  multiplicative  uncertainty  in  Figure  5.4.  The  uncertainty 
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weighting  function  Wn  is  assumed  as 


Wn  =  0.0045 


s/0.09+1 
s/80  +  1 


(4.15) 


The  magnitude  of  the  uncertainty  block  |A|  is  bounded  by  one.  The  ellipsoids  to  contain  the  reachable 
set  are  calculated  based  on  quasi-LPV  models  {Q-LPVj  and  Q-LPVf)  and  shown  in  Figure  5.5.  It  is 
noticed  that  the  function-substitution  quasi-LPV  model  has  smaller  ellipsoid  to  represent  the  reachable  set 
than  the  Jacobian-linearization  quasi-LPV  model.  Recall  that  the  state  X2  is  measurement  for  design  an 
output  feedback  LPV  controller.  It  is  noticeable  that  the  axis  of  the  ellipsoid  at  X2  direction  is  smaller 
of  the  function-substitute  quasi-LPV  model  than  the  Jacobian-linearization  quasi-LPV  model.  The  set  of 
{(xi,  X2)|x^Px  <  1}  in  Figure  5.5  is  much  larger  than  the  set  shown  in  Figure  5.3  for  each  quasi-LPV  model, 
respectively.  It  is  obvious  that  adding  unmodel  dynamics  in  the  quasi-LPV  model  enlarges  the  size  of  the 
set  {(xi,  X2)|x^Px  <  1}.  The  state  xs  of  augmented  open-loop  system  is  the  state  of  the  weighting  function. 
The  sizes  of  the  ellipsoids  in  X3  direction  are  similar  to  each  other,  based  on  each  quasi-LPV  model  {Q-LPVj 
and  Q-LPVf). 

In  this  example,  adding  unmodel  dynamics  in  quasi-LPV  models  does  not  change  the  comparison  result 
that  the  function-substitution  quasi-LPV  model  is  less  conservative  to  present  the  nonlinear  dynamics  than 
the  Jacobian-linearization  quasi-LPV  model.  Note  that  it  is  unknown  that  adding  an  uncertainty  block 
changes  the  comparison  results  in  general. 


4.3.2.  Parameter-Dependent  Lyapunov  Function.  A  Lyapunov  function  V (x)  =  x^P(x)x  is  used 
to  calculate  the  smallest  set  {x|V(x)  <  1}  which  contains  the  reachable  set  of  the  nonlinear  dynamics. 
Consider  quasi-LPV  models  without  an  uncertainty  block.  To  solve  the  LMI  optimization  problem  in  equa¬ 
tion  (3.17),  basis  functions  of  P(x)  are  required.  A  matrix  function  P(x)  and  the  time  derivative  of  P(x) 
can  be  written  as: 

N 

P{x)  =  Pi  e  P""""  (4-16) 

i=l 

=  (4.17) 

j=l  i=l  ^ 


with  given  basis  functions  fi  (x) ,  where  s  and  N  are  number  of  scheduling  parameter  and  number  of  basis 
functions,  respectively.  In  the  example,  the  set  of  basis  function  {/i(x)}  is  defined  as  the  first  order  polynomial 
function  set  {1,X2}  for  computational  convenience.  The  time  derivative  of  P(x)  is  written  as: 


P(x)  =  X2P2 


(4.18) 


where  P2  G  ,  In  conventional  parameter-dependent  Lyapunov  function  LPV  synthesis  [3],  constant 

bounded  values  v  of  the  parameter  rates  are  used  to  describe  the  time  derivative  of  P{x)  as: 


N 


P{x] 


=  E(E±-. 

j=i  i=i 


^  dXi 


<  jyj 


(4.19) 


where  ±iyj  represents  all  possible  combination  set  of  —lyj  and  jyj.  In  this  paper,  the  time  derivative  of  P(x) 
is  written  as 


n  N 

p(^)  = 

j=l  i=i 


dfijx) 

dxj 


(4.20) 
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where  gj{x)  :  IZ  is  satisfied  with  \xj\  <  gj{x).  To  evaluate  equation  (4.20),  the  function  gj{x)  is 

estimated  from  nonlinear  dynamic  simulation  results  with  the  bounded  energy  inputs.  For  this  example, 

g{x2)  =  (1  +  f{x2))x2  (4.21) 


where  f{x2)  is  defined  in  equation  (4.6). 

The  boundary  of  the  set  {x|x^P(x)x  <  1}  to  contain  the  reachable  set  IZ^e  is  shown  in  Figure  5.6.  The 
solid  line  in  Figure  5.6  represents  the  boundary  of  the  set  {x|x^P(x)x  <  1}  calculated  based  on  the  function- 
substitution  quasi-LPV  model.  The  dotted  line  in  Figure  5.6  represents  the  ellipsoid  calculated  using  the 
singular  quadratic  Lyapunov  function  based  on  the  function-substitution  quasi-LPV  model.  It  is  noted 
that  the  set  {x\V{x)  <  1}  using  the  parameter-dependent  Lyapunov  function  is  smaller  than  using  singular 
quadratic  Lyapunov  function.  This  result  corresponds  to  that  LPV  control  synthesis  methodology  using  the 
parameter-dependent  Lyapunov  function  leads  a  less  conservative  result  than  using  the  singular  quadratic 
Lyapunov  function  [3].  The  dashed-dot  line  in  Figure  5.6  represents  the  set  {x|x^P(x)x  <  1}  calculated 
based  on  the  Jacobian-linearization  quasi-LPV  model.  It  is  observed  that  the  function-substitution  quasi- 
LPV  model  is  less  conservative  than  the  Jacobian-linearization  quasi-LPV  model  to  present  the  reachable 
set  of  nonlinear  system.  The  analysis  of  quasi-LPV  models  with  an  uncertainty  block  leads  the  same  result 
that  the  function-substitution  quasi-LPV  model  is  less  conservative.  The  plots  of  the  calculated  sets  are 
omitted  for  space  limitation. 


4.4.  LPV  Controller  Synthesis.  An  LPV  controller  is  synthesized  based  on  each  quasi-LPV  model 
and  is  simulated  with  the  nonlinear  simulation  to  compare  the  closed-loop  system  performance. 

The  prime  objective  of  an  LPV  control  synthesis  is  to  track  a  given  command.  An  LPV  controller 
synthesis  formulation  is  taken  from  Ref.  [2]  shown  in  Figure  5.7.  The  performance  weighting  function  Wi 
and  control  weighting  function  W2  are  taken  from  Ref.  [2]  as 


Wi 


0.5 

s-h  0.002’ 


W2 


0.1s 

s  -h  1000 


(4.22) 


The  noise  weight  function  is  defined  as  constant  0.05  over  all  frequency  range  to  present  5%  measurement 
error. 

To  solve  the  LMI  optimization  of  LPV  controller  synthesis,  a  scheduling  parameter  is  discretized  over 
all  spaces.  The  grid  points  are  presented  in  equation  (4.13).  The  LPV  controllers  are  synthesized  for  each 
quasi-LPV  model  with  same  weight  functions,  using  standard  software  from  the  MATLAB  LMI  toolbox  [12]. 
There  are  two  generated  LPV  controllers  (Q-LPVj,  Q-LPVf)  according  to  two  quasi-LPV  models. 

The  responses  to  step  change  in  command  from  -1  to  5  of  the  closed-loop  system  are  shown  in  Figure  5.8. 
It  is  observable  that  the  overshooting  responses  of  the  LPV  controller  synthesized  based  on  the  function 
substitution  quasi-LPV  model  are  smaller  than  those  with  the  Jacobian-linearization  quasi-LPV  model.  The 
calculated  7  in  equation  (2.3)  is  0.92  for  the  Jacobian  linearization  quasi-LPV  model  and  0.73  for  the  function 
substitution  quasi-LPV  model. 

The  performance  results  are  correspond  to  the  analysis  results  of  LPV  models.  The  ellipsoid  based  on 
the  functional  substitution  quasi-LPV  model  to  present  a  reachable  set  of  the  nonlinear  system  is  smaller. 
That  implies  that  the  quasi-LPV  model  generated  by  the  function  substitution  is  less  conservative.  Note 
that  it  is  not  known  that  the  function  substitution  method  can  always  generate  a  less  conservative  quasi-LPV 
model  in  terms  of  an  ellipsoid  which  contains  the  reachable  set. 
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5.  Conclusion.  In  this  paper,  one  of  approaches  to  compare  quasi-LPV  models  which  represent  a 
nonlinear  system  is  demonstrated  in  terms  of  the  smallest  set  which  contains  the  reachable  set  of  the  nonlinear 
system.  Based  on  the  size  of  the  set,  it  is  possible  to  define  which  quasi-LPV  model  is  less  conservative  to 
present  the  reachable  set.  The  quasi-LPV  models  of  a  nonlinear  system  are  generated  by  three  different 
methods  (Jacobian  linearization,  state  transformation,  and  function  substitution)  to  facilitate  to  design  an 
LPV  controller  of  the  nonlinear  system.  LPV  controllers  are  synthesized  based  on  each  quasi-LPV  model 
and  simulated  with  the  nonlinear  system  to  compare  the  closed-loop  performance.  The  performance  results 
correspond  to  the  results  of  comparing  the  size  of  the  set  calculated  based  on  each  quasi-LPV  model.  Based 
on  the  analysis  results  of  quasi-LPV  models,  it  is  possible  to  choose  which  LPV  model  is  used  for  LPV 
controller  synthesis  of  a  nonlinear  system. 
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Input  case  1 


time 


Fig.  5.1.  Input  signals. 
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